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ABSTRACT 


A geometrically similar family of three supercavitating 
hydrofoils was tested in the M.I.T. variable pressure water 
tunnel. The half-span foils were of elliptical planform; 
the ratio of foil half-span to tunnel height was 1/4 for the 
Small foil, 1/2 for the medium foil, and 3/4 for the 
largest foil. The tunnel was of square cross section. Lift, 
drag, moment, tunnel speed, ambient pressure, and cavity 
pressure were measured for attack angles from 8 to 21 degrees 
and a variety of ambient pressure settings; cavity length 
measurements were obtained from photographs. Results were 
compared with theoretical results obtained by Leehey, and 
also with a more detailed numerical lifting surface procedure 
developed by Jiang and Leehey. For the small and medium foils, 
it was sufficient to correct only for the effect on downwash 
of the images of the trailing vortices. The large foil data, 
however, required further correction; upon application of the 
same corrections which were applied to the data for the two 
smaller foils, the force and moment data for the large foil 
plotted slightly lower than did the data for the two smaller 
foils, while the cavity length data for the large foil 
indicated cavity lengths significantly larger than for either 
of the theoretical predictions or the cavity length data for 
the two smaller foils. Through the application of existing 
two-dimensional corrections, the force data for the large 
foil were brought into close agreement with the force data 
for the two smaller foils, but no suitable correction factors 
exist for the cavity length data. 


Thesis Supervisor: Patrick Leehey, Ph.D. 
Title: Professor of Naval Architecture and 
Applied Mechanics 


TOAATIS SILK TO GTS 
TASS APiv14 So Sr lovee 


9,18. Tei leah eae 

ZWolel wolves ele WS 
4 (S$2a83 Solan (7% ‘ a 
near. | 108 CGR ‘ 


len eae 


£3 
ACKNOWLEDGMENTS 


This research was carried out under the Naval Sea 
Systems Command, General Hydromechanics Research Program, 
Subproject SR 009 01 01, administered by the David W. Taylor 
Naval Ship Research and Development Center, Contract 
Number N00014-76-C-0358. The author is exceedingly grateful 
to Professor Patrick Leehey of the Department of Ocean 
Engineering at the Massachusetts Institute of Technology 
for his inspiration, suggestions, and guidance. Special 
thanks are also due to S. Dean Lewis and C. W. Jiang, also 
from M.I.T., who rendered valuable assistance in the 
performance of the experiments, and to Miss Cheryl Gibson 
for typing the various drafts and for enduring countless 
errors. Finally, the author is indebted to his wife, 
Gretchen, for her unfailing support and encouragement 
and for the many sacrifices endured by her during the 


preparation of this report. 


~ 


i -6e [ere £9 15 °00 £8 


ee ae eaw : 
el nerle morta Jf3sre@ 
’ = 
uqOlsve@ Bas sotesagie 
oH? NEE0-9-ay—p ee 

a retired ¥ b3ec9 
J-830733s60! oid ga 
sees 5 . oottesigual ae 
5 ols gaa 


TABLE OF CONTENTS 


IRE PAIGE Misislcialscleiscc ccs ect ce ste cessesce Sretalcionerene oielele: oie 
ene PME MOMMAY oie7e''o)/s/561c 610%! ©) 6, 0) «6 © 1c we) ofl oe! ai 6-6: are % 8 ol oo e8 1 gus oe 
EMIMINETIOGMEN TS «so 00s occ ec esc eeeeseeseeesvccccccessee 
TABLE OF CONTENTS ....... Sc) se loons 5 oO CEE ee re 
foe OE FIGURES ....cccecee sholetoanetatetel seltstelis| 6iels vets Sy ereis ce) « 
NOMENCLATURE ...... gua lee: ole yel craton si chonsrenet forbs 6s oS. wo eee wiy wwe 
PME RODUCTION .cccccccccscesccocsece EWeiclictes clio ts sicellea\@) sens 
EXPERIMENTAL APPARATUS AND PROCEDURE ......2.cccccsees 
APPAFACUS Sese ccc wes mets ese sece svc sccencs cs eres 
AO COC Utse mearaneter el clciicte Giisl's!/ie\ 0 clo © 6 01 e's! s ei oie 0 See's 6 0s0 8 
POPS GIS) 4 Giga > C016 516 Occ O1O OIG CROICICIORO IOC I CHORE Ra CaCHr arin ar armrests ee 
Pein: so ciersuenaieis) oie S18 GCs OIG GORE ONCE IDLO ODE © Oa ete siels 
DEAG c's « cleisic aie 6 06 oso ce cee ee alesis ole 616616 we 6 6 6's Sistas 
MOMENE .ecccsccseccceccscce Srchialieeirelaiiel sie ole ‘oie0 Siie (aes 
CAVAUCY GENGEN csc cesseccss aati au suel eters sPotalaletecoiers 


CAVEEACHONONUMDEGL occ ccccsseccavesccnetiseces Sp Soscekee 


DES CUSSION «6.2 ccc cc eee cocccns ess eee we 0) eiei6 «70 SOS areca 
PONCLUSIONS 5. ccc wees ccccccccccseses Cie siei's) stele (se oie. es wie 
BREPERENCES ........ sistlelielistielis) o/s siere «e's es s\'e)s\ © 210s ee sie) «)'6\0 6 610 


APPENDIX I ... eeesveee ceoeeeeaeeeeeeeeveeeeeeeeeteeeeeeeeee eeeve 
APPENDIX agg aeee50xeee8ee eeeeveeeee eoeeeeoeeeeeeeeeeee eeeeeee 


PUI LG oc ce cc cee pete cece sigs) elses! © © sisies) 016 66 5106 
APPENDIX IV @eeeesoesvcseseerteeeeeeeeeee eeeec0seeeeeeeeeeteeeeeee @ 


Va 


** @€s2)G8 464% + <4 obat@e 


‘= @ we « ‘ 
« ® ° « 
. 
J 

* * 

} 
> VPlpy 

=a) 


$8432 10 sonbgd 


-#*+ en €64@ 


""-'e-2 as Oe ae 
a 


- 


-+e= © & a ea 
, 

FY ‘ — - ee ee ee . a if 
eae < an 

27 ee ee ee : 


"See eae oe = 


TPR Meee Its) sap eas Swag) . 


Sree 
$204 /us 6S graced 


— 


Figure 


a 
2 
3 


LIST OF FIGURES 


Mor. t. Test Faciirty 


Foil Dimensions 


Cy / Om vs. CO / Ons 


Cy/ Om vs. OT /One 


eeeee 


eocoeoeeeveeseeeeeeeeeee @& 


a = 80% 
a= QPS? 
a dU leo 
a ZO 
a Ny aoe 
a LG 20° 
Cie — eke 207 
a 21 20° 
P= 48.0 
G. =79..5% 
a 0 
a 12: 0'° 
a 14.0° 
a 16.0° 
a 30° 
a 21.0'> 
a 8.0° 
q@= 9.5° 
1 a0 


coeoeoeveeeeeeeeeeeee ee @ 


oeeeeeveeeeeeeeeeeeee 


eoeeeeeeeeeeeeeeeee @ 


ceeceeeeeeeeeeeee eee 


eoseeeeeeeeeeeeeee eee 


eoeeeoeeeeeeeeeeeeeee 


eeeeoeneneeeveeeeeeeeeeee 


ov Gekiiost tee? 
eta te ena/ aoomeg 
oo TM 2xgss98F tik 


Gad 


ui imgd\,.2 «ev 


~ 


Figure 

23 Cy/ On vs. O/ ome CELA SO 5 ieee ee ae 
24 Ce eA Ob etait eis Ges aio Ga dw o8 <4 
25 Gide = Pale ONION etre encoh eye iene. vrciss..¢. 0.l000's, ete 6 
26 (ST 118} (0) lee a eae a 
27 Cee ele Oe eahie le,°ey 0 '4:-ove (0.78 ee! ei ete wus 
28 0 / On NR mle a seh O Ils oe qucceuels seis: 6, uo 69 6 6 6% is 
Zo MOG E Ol pei cnetensc.0.6 6 0.5.0 di 6.0 « o 6 
30 rae Ge Ph Oa ISS Aid eet a shes. eis a etek es) 
onl Comparrson or "Cavrty “Lengens 25 shoes oo. SF Meee Ss 
32 O, VS. OL, SMe Oa tarctats sratenete ce eres spepens co e.ehe 6s 
33 MC CMsM er Odal- sisle.s is abave. eye ienelt6 6.6 sensei ec 
34 Se cto Pay 2c )sh DEER Re ea eee ae loi bierohnio aes 
35 Large Foil: Comparison of Cavity 


Pressure Measurement 
TEechnagques. .<...6.e sss Ges 


& (2) 


NOMENCLATURE 


fou Leaaing edge bevel Tength; see Fig. 2 
aspect ratio = 2s*/s 

foil mean chord, measured at centroid position 
foi) root chord 


tunnel correction factor for influence of images 
of trailing vortex system 


drag coefficient = t0ta’ drag = c, + ACD 
(1/2) pu2s uncor i 


uncorrected drag coefficient, as calculated from 


the raw experimental data = measured rag 
(1/2) pU“S 


drag coefficient, corrected for blockage; see 
equation (1) 


change in induced drag due to effect of images 
of trailing vortex system 


measured lift 


lift coefficient = (1/2) pu2s 


two-dimensional lift coefficient; see equation 
(II-5) 


moment coefficient about the midchord 


point = moment 
(1/2) pU2Sc 
tunnel height = 51 cm 


cavity length as a function of spanwise position; 
measured from leading edge 


cavity length at centroid position; measured 
from midchord 


minimum pressure along tunnel wall in vicinity 
of cavity 
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Measured cavity pressure 


upstream tunnel static pressure at the spanwise 
position of the foil centroid 


vapor pressure of water 

tunnel static pressure controller setting 
static pressure signal to controller from 
upstream, downstream tunnel static pressure 
taps; see Fig. 3 


strength of discrete source located in the 
(i,j)-th element 


source strength at (x,z) 

foil half-span 

element length for (1i,j)-th element 

planform area of half-span model 

tunnel test section cross-sectional area 
thickness of hydrofoil 

perturbation velocities in the x,y-directions 


water velocity one meter upstream of dynamometer 
shaft axis 


maximum water velocity (at point of minimum 
pressure, Pp) ievicanity of cavity 


chordwise coordinate of leading edge as a 
function of spanwise position 


distance of cavity pressure tap from midchord 
foil coordinates; see Fig. 2 


distance of centroid of half-span foil model 
from midspan 


distance of cavity pressure tap from midspan 
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y (x,2) 


geometric angle of attack with reference to 
tunnel side wall; when used in Appendices II 
and III, a indicates angle of attack in a free 
stream 


"true" angle of attack, as corrected for effects 
of images of trailing vortices = a + Aa; 


change in angle of attack due to effect of 
images of trailing vortices 


arctan (2a/o,) 


strength of discrete vortex located in (i,j)-th 
element 


vortex strength at (x,z) 


Tatsao of. fol frontal width. to, tunnel, width; 
see equation (3) 


ratio of foil frontal area to tunnel cross- 
sectional area (without foil); see equation (4) 


dummy variables in the x,z-directions 
density of water 


cavitation number based on measured cavity. 


pressure = (1/2) pu? 


Cavitation number, based on vapor 


(1/2) cuUZ 


cavitation number, corrected for blockage; see 
equations (2) and (5) 


cavitation number, based on Pp: P and V, 


rel 
the maximum water velocity = —Pb™Pc_ 
(1/2) pv2 


od soueiniss siscw sone ae 
Tl ecolionggA al beer ASAY 41 
—> aie | 8 =A45 3 as3 Sita yeas *% 


i P 


o iit 


t kee 
on 


C LY. pam af 
{ 
alo Ie 
‘ 

6 
[y 

! 

be —_ 
7 


aaa 
oJ f.0 TON. eee > a 


> 


a= 


INTRODUCTION 


Corrections for water tunnel wall effects during the 
testing of subcavitating hydrofoils are essentially the 
same corrections developed fifty years ago for use in wind 
tunnel testing of airfoils; references [1] and [2] are 
standard texts in this area. Only within the last twenty 
years have wall effects been considered in the case of 
supercavitating hydrofoils, with the majority of theoretical 
and experimental work concentrated on two-dimensional hydro- 
foils. The problem of blockage can be significant; the 
shape of the foil and cavity combination must be known in 
order to apply a blockage correction of the type utilized in 
subcavitating flow. Unfortunately, this proves to be an area 
-Of great difficulty in the various theories. Corrections to 
drag coefficient and cavitation number for the two-dimensional 
pure-drag case are given in references [3] through [7]; a 
nonlinear theoretical model of the two-dimensional lifting 
problem is formulated in reference [6]. In reference [8], 
Baker converts this model into a computer program which 
numerically calculates the wall effect on two-dimensional 
hydrofoils of arbitrary shape. It was hoped by Baker that 
these results would prove adequate for application to flows 
Over high aspect ratio hydrofoils. Prior to the current 


research, a systematic series of experiments which exhibit 
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wall effects on a geometrically similar family of hydrofoils 
had not been carried out in either the two- or three- 
dimensional case. The purpose of the present research is 

to make available data on a geometrically similar family of 
three supercavitating hydrofoils of finite span, and to apply 
existing two-dimensional corrections in an attempt to bring 
the data into agreement with three-dimensional, unbounded 


flow theory. 
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EXPERIMENTAL APPARATUS AND PROCEDURE 


Apparatus 


The experiments were conducted in the M.I.T. recircula- 
ting, variable pressure water tunnel (Fig. 1). The test 
section is 51 cm square in cross-section, and 147 cm long. 
Test section velocity was continuously variable between 0 
and 9.1 m/sec; a manometer indicated the difference in 
pressure across a contraction in the tunnel one meter up- 
stream of the dynamometer axis. This manometer had been 
previously calibrated using a pitot tube in the test section. 
The assumption was made that the blockage caused by the foil 
and cavity combination did not affect the manometer 
appreciably. Test section static pressure was variable 
between 60 mm Hg absolute and atmospheric pressure; static 
pressure readings were taken using a mercury manometer 
connected to either or both of two taps on the centerline 
of a side wall of the test section. These taps were 
located 62 cm upstream and downstream of the hydrofoil 
dynamometer, which was located at the center of the test 
section. 

Foil dimensions are given in Fig. 2. The three half- 
Span, aspect ratio five foils had sharp leading edges. 
Geometrically similar elliptical planforms were chosen even 
though foils of this type were more difficult to fabricate 


than, for example, foilsof rectangular planform; elliptical 
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planform foils exhibit less tip loading effect than do 
rectangular planform foils, thereby rendering them more 
amenable to analysis. Each foil was drilled with a hole 
terminating in a pressure tap on the cavity side of the 
foil. The medium-sized foil was the aspect ratio five foil 
utilized by Leehey and Stellinger [9], except that the 
cavity pressure tap was relocated to a position nearer the 
tip. The composition of the medium foil was of type 304 
Stainless steel, while that of the small and large foils 
was of type 614 aluminum bronze. Each foil was welded to 
its own 10.16 cm diameter mounting plate; the mounting 
plates, in turn, were bolted to the dynamometer so that 
the small, medium, and large foils protruded 1/4, 1/2, and 
3/4 of the way, respectively, downward from the top window 
into the test section. 

The dynamometer was capable of measuring three 
moment and three force components, but only the moment 
about midchord and the forces tangential to and perpendicular 
to the span and chord were measured. Zero geometrical 
angle of attack was set by aligning the flat (wetted) side 
of each foil with the test section side wall. Upon 
alignment, a sighting telescope, which was attached to the 
dynamometer, was aligned with a scale on the laboratory 
wall; this scale, graduated in increments of 0.05 degree, 


was utilized in setting the geometrical angle of attack. 
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Procedure 
Before test runs were begun on any one day, the tunnel 
was operated at medium velocity for twenty minutes at a 
pressure of 300 mm Hg absolute; this removed almost all of 
the air trapped in remote locations of the tunnel. Prior 
to each test run, the bubble chamber shown in Fig. 3 was 
utilized to collect and draw off any gases or air which may 
have been trapped in the static pressure sensing line; once 
Cleared of gases, the bubble chamber was merely a reservoir 
of common pressure. paeene other hand, small amounts of air 
were bled intermittently into the cavity pressure sensing 
line during testing in order to keep it free of moisture. 
Each foil was tested at a variety of attack angles 
from 8 to 21 degrees. The lower limit of 8 degrees was 
chosen so that the cavities would not terminate on the foil. 
For all test runs, the test section velocity was set at 
8.5 m/sec, whereupon static pressure was reduced until a 
cavity covered the entire planform. Readings were taken 
of lift, drag, moment, static pressure (with pressure taps 
U and D open; see Fig. 3), cavity pressure, and velocity. 
A Polaroid photograph was taken of the foil and cavity from 
the wetted side of the foil, utilizing an exposure time 
of 1/40 second with bottom flood lighting. Pressure tap D 
was then closed, and the same readingstaken; reopening 
pressure tap D and closing pressure tap U, the same series 


of readings was repeated once more. (Photographs were not 
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taken for these last two pressure tap combinations since 
hardly any variation in cavity shape or length was noted 
when compared to the situation where both pressure taps 
were open.) This completed the measurement of three data 
points at one controller setting. Both pressure taps 

were then opened, and the static pressure was reduced by 
20 to 30 mm Hg, whereupon the above procedure was repeated. 
The static pressure was reduced in this fashion until an 
excessive number of recirculated cavitation bubbles resulted 
in either extremely poor conditions for photography or 
rapidly fluctuating instrument readings. Prior to and 
immediately following a run of this type, the room and 
water temperatures were recorded, as were the tares on the 
instrumentation. 

A computer program was utilized in the data reduction. 
Wetted surface frictional drag (for each foil and its 
mounting plate) and dynamometer shaft twist were taken into 
account, and test section static pressure readings were 
modified for static head difference to give the static 
pressure at the vertical position of the foil centroid. 
(Gravity effects were otherwise neglected.) Cavitation 
number was computed using measured cavity pressure, and 
also by using vapor pressure at the ambient water temper- 
ature (computed using equation (2) on page 151 of 


reference [10]). Lift, drag, and moment coefficients 
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were also calculated. Standard wind tunnel corrections 
(as described in Appendix I) were applied to the data to 
account for the effect of the images of the trailing 
vortices; this resulted in an increase in drag coefficient 
and also in an increase in effective angle of attack. 
Utilizing the photographs, cavity lengths were measured 


from the midchord at the centroid position. 
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RESULTS 


Experimental results were compared with predictions 
from Leehey's linearized theory [11], which utilizes the 
method of matched asymptotic expansions; the theory is valid 
to first order in angle of attack, and to second order in 
the reciprocal of aspect ratio. Comparisons were also 
made with a more detailed linear lifting surface calculation 
developed by Jiang and Leehey [12]. This numerical theory 
is expected to be more accurate at the lower values of 
0 /Ons where the cavity length is the same size as the foil 
span or longer; it should also be more accurate, in general, 
for prediction of moment coefficients. Synopses of the 
theories developed in references [11] and [12] are given 
in Appendices II and III, respectively; the experimental 
data are given in Appendix IV, both in raw form and with 
standard wind tunnel corrections (Appendix I) applied. 

The wall boundary layer momentum thickness was 
approximately 1.9 mm, as inferred from measurements by 
Stellinger under similar conditions [9]. Since less than 
two percent of the wetted surface area of the small foil 
was within 1.9 mm of the upper wall, the upper wall 
boundary layer was assumed to have negligible effect on the 
force and moment coefficients of all three of the foils 


tested. 
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It should be noted that the only correction for wall 
effect which has been applied to the data in Figs. 4 
through 30 is the correction for the effects of the images 
of the trailing vortices (Appendix I). Also worthy of note 
is the fact that the supercavitating condition was never 
achieved for the small foil at an angle of attack of 8 
degrees. This fact bears upon the lack of scaling evident 
in the force and moment data of Figs. 4, 12, and 20 for this 
foil. This matter will be considered again in the sub- 


sequent discussion section. 


Lift (Figs. 4 through 11) 


There was very good agreement, overall, between 
theoretical predictions and the experimental data. For 
small values of the similarity parameter F/O (7eu,+ £or 
long cavities), there was much better agreement with the 
numerical theory than with the asymptotic theory; in this 
case, the foil and cavity combination no longer has a 
large aspect ratio. For the higher values of O /Ons 
Leehey's theory appears to underpredict, especially at the 
higher angles of attack. This is in apparent contradiction 
with Fig. 4 of reference [9]; it should, however, be noted 
that the data of reference [9] were based on cavitation 
numbers calculated with vapor pressure rather than with 


measured cavity pressure. It can be seen from Figs. 32 
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through 35 of the present work that, when based on measured 
cavity pressure, the cavitation number will be smaller, so 
that the above contradiction is only apparent. The small 
and medium foil data plotted almost on top of one another, 
whereas the large foil data plotted lower than did the small 
and medium foil data, especially at attack angles greater 


than 11 degrees. 


Drag (Figs. 12 through 19) 

When the effects of streamwise foil or flow curvature 
are negligible, linear theory predicts Cy = OmnCy so that 
the drag data are plotted as Ch/an vs. 0 o/s Ch/an plotted 
somewhat higher than the theoretical predictions except 
for very long or for very short cavities. As with the lift 
data, the small and medium foil data plotted together, with 
the large foil data plotting somewhat below these, the 


effect becoming more pronounced at the higher angles of 


attack. 


Moment (Figs. 20 through 27) 


The moment coefficient was taken about the midchord 
consistent with the right-hand rule. As pointed out by 
Leehey and Stellinger [9], Leehey's matched asymptotic 
expansion theory neglects lifting surface corrections 
(third order in the reciprocal of aspect ratio) for the 


moment coefficient. It is not surprising, therefore, that 
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the current experimental data agreed best with the more 
detailed lifting surface theory of Jiang and Leehey [12]. 
The load cell utilized for moment measurements on the small 
foil was of inadequate sensitivity, and consequently, much 
of the small foil data at the smaller attack angles may be 
less accurate than the data for the medium and large foils; 
comparisons must therefore be made between the large and the 
medium foil data. Once again, the large foil data plotted 
below that of the medium foil, especially at the larger 


attack angles. 


Cavity Length (Figs. 28 through 31) 


As mentioned previously, cavity lengths were measured 
from midchord at the centroid position rather than from the 
leading edge, since this convention was utilized by Leehey 
L111]. Cavity length (L) was nondimensionalized on the foil 
mean chord (c). The cavity length data for the small and 
medium foils plotted on top of one another (Figs. 28 and 29); 
in both cases, the cavity lengths were slightly less than 
those predicted by theory, but the overall agreement with 
theory was good. In contrast, there was a marked deviation 
in the large foil cavity length (Fig. 30) and shape (Fig. 31), 
especially at the root, where local blockage was the worst. 
The tip vortices seen in the photographs had no evident 


effect upon the force or moment coefficients. Neither 
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Leehey's matched asymptotic expansion theory nor the numer- 
ical lifting surface theory developed by Jiang and Leehey 
accounts for the roll-up effects which produced these 
isolated trailing vortices in the wake. 

Representative photographs of the large, medium, and 
small foils at the same angle of attack and at approximately 
the same cavitation number are shown in Fig. 31. The cavity 
shapes shown were averaged over the exposure time of 1/40 
second; visual observations showed that the instantaneous 
cavity shape was highly unsteady. For the small and medium 
foils, the general outline of the cavity was elliptical, 
as predicted by theory; the cavity length was undoubtedly 
reduced toward the tip due to the hydrostatic pressure 


gradient existing in the tunnel. 


Cavitation Number (Figs. 32 through 35) 


Below a cavitation number of approximately 0.4, there 
was almost no difference between cavitation number calculated 
with measured cavity pressure or with vapor pressure. It 
should be noted that in this region there exists a large 
variation in cavity length from very short to very long. 

This indicates that for even short cavities, where ram 
effects of the re-entrant jet from the cavity end have 
heretofore been thought to have a pronounced effect on 
measured cavity pressure, the actual cavity pressure was 


approximately equal to the vapor pressure of the water. 
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Above a cavitation number of 0.4, the disparity between 
Oy and G.. increased with increasing cavitation number; the 
data in this region, however, were only for large angles of 
attack and small cavity lengths, indicating that there may 
indeed have been some effect caused by the impinging of the. 
re-entrant jet on the cavity pressure tap. 

After the series of basic experiments was completed, 
ram effects were further investigated on the large foil 
by Jiang and Leehey, this time utilizing a total head tube 
for cavity pressure measurement. The L-shaped total head 
tube protruded downwards into the cavity from the upper 
tunnel wall so that it was parallel to the foil wetted 
surface and pointed towards the leading edge and away from 
the impinging re-entrant jet. Cavity pressure readings 
were taken with both the standard foil surface pressure tap 
previously utilized and the total head tube, and vapor 
pressure was calculated as was done previously; results are 
shown in Fig. 35. For the larger attack angles and shorter 
cavities (higher cavitation numbers), use of the total head 
tube significantly reduced re-entrant jet effects; the 
result was that cavitation numbers calculated with total head 
tube data were as much as 5 to 8 percent higher than 
cavitation numbers calculated with data obtained from the 
previously used foil surface tap. Although a similar check 
was not made on the small and medium foils, it can be 


assumed that the same general trend would occur. It should 
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be noted, however, that data points in Figs. 32 and 33 
exhibit a more pronounced downward "hook" at higher 
cavitation numbers than do the large foil data (Figs. 34 
and 35), a fact which accounts for the noticeable "hooks" 
in the lift and moment data at the higher angles of 

attack (see, for example, Figs. 10, 11, 26, and 27). 

If the total head tube had been utilized throughout the 
entire series of experiments, these deviations from theory 
at high attack angles and short cavities would probably 


have been less. 


Influence of Static Pressure Controller 

As stated earlier, readings were taken for the various 
combinations of the two static pressure taps; this was 
done in order to determine the effect of blockage on 
pressure and velocity at the two different tap positions. 
It was anticipated that this variation in static pressure 
tap combination would have no effect on the forces and 
moments experienced by the foil. Due to the location within 
the system of the static pressure controller (Fig. 3) 
(which was inadvertently overlooked), the forces and moments 
varied significantly when the combination of static pressure 
taps was changed. The explanation for this is as follows. 
Before data were taken, both static pressure taps were opened 
and the controller setting reduced to p; the actual static 


pressure in the tunnel was decreased by the control system 
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until the combined signal from both taps was p. When the 
instrument readings settled out, data were taken. The 
downstream static pressure tap was then closed, so that the 
only signal to the controller was Py: which was greater than 
the controller setting, p. The controller therefore lowered 
the actual static pressure in the tunnel until Py equaled p. 
Consequently, the effective cavitation number (as seen by the 
foil) decreased, causing the cavity to grow. This reduced 
the effective camber of the foil and cavity combination; 
this loss of effective camber resulted in a reduction of the 
force and moment coefficients from their previous values. 
On the other hand, the measured cavitation number increased 
Slightly due to a small decrease in measured cavity pressure; 
the presence of the controller prevented the measured static 
pressure from changing significantly. The opposite effect 
was observed with the upstream tap open and the downstream 
tap closed (i.e., an increase in force and moment coefficients, 
an increase in effective cavitation number, and a slight 
decrease in measured cavitation number when compared to the 
case when both pressure taps were open); there were only a 
few of these data points taken. 

Although the desired pressure readings were not obtained 
in the current experiment (i.e., the actual tunnel static 
pressure was Gifferent for each of the tap combinations), 


the results do show the significance of the blockage. When 
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comparing the results obtained with the different tap 
combinations for the small and medium foils, there was very 
little change in the plotted data points; the increased 
blockage in the case of the large foil resulted in signi- 
ficantly different force and moment data. For the sake 

of clarity, the data for the different tap combinations are 
shown only for the large foil at an angle of attack of 

14 degrees. Unless otherwise specified, further remarks 
concerning foil data will be for the data obtained with both 
static pressure taps open. 

The general trend of the data agrees with the experi- 
mental results obtained by Kermeen (See page 37 of reference 
[13]); he observed that the longer the cavity in super- 
cavitating flow, the smaller the force coefficients. The 
same conclusion follows from the theoretical results of 
Leehey, as may be seen from Fig. 4 of reference [9]. 

At a later time, the system indicated by the dashed 
lines in Fig. 3 will be utilized to obtain data on the large 
foil; the present system would be operated with only tap U 
open to provide an indication of static pressure at "infinity" 
to the controller. This configuration will allow pressure 
measurements to be made without affecting the input signal 


to the static pressure controller. 
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DISCUSSION 


For the small and medium foils, it appears to be 
necessary to correct only for the effects of the images of 
the trailing vortices in accordance with standard wind 
tunnel procedure. These corrections to drag coefficient 
and effective angle of attack brought the data into good 
agreement with theory; the lift, moment, and cavity length 
data also agreed well with theoretical predictions. For 
the most part, the plots for the small and medium foils 
were very close to each other. 

Upon application of the same corrections to the large 
foil data, good agreement with theory was seen for the force 
and moment coefficients; they were, however, generally 
lower than the data for the smaller foils, especially at 
higher attack angles. The plot of cavity length for the 
large foil, however, showed substantially longer cavities 
than predicted by either of the two theoretical models used 
for comparison. The most pronounced deviations in cavity 
length occurred for the long cavity (lower oon) data. 

When compared with the data obtained from the two 
smaller foils, the general trend of the large foil data is 
in agreement with Baker's predictions for wall effects on 
two-dimensional supercavitating hydrofoils [8]. Baker makes 


the point that the wall effect is negligible for tunnel 
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height-to-foil chord ratios (H/c) greater than 10 for 
cambered foils. Although he makes no such generalizations 
concerning supercavitating flat plates, (i.e., the case 
applicable to the current research), it is seen from 

Figs. 4 through 27 that the effect of blockage on the force 
and moment coefficients is not exceedingly large. The 
following line of reasoning explains why the force and 
moment data for the large foil were somewhat lower than for 


the two smaller foils. For the three foils at the same 


cavitation number and at the same angle of attack, the 


Magnitude of the velocity on the cavity wall of all three 


foils is 0/2 (from linearized theory); consequently, all 
three foils have identical pressure coefficients on the 
cavity. The largest foil experiences significantly more 
blockage than does either of the smaller foils, so that the 
velocity on the wetted surface of the large foil is much 
higher than for either of the smaller foils; the pressure 
coefficient on the wetted surface of the largest foil is 
therefore the least of all three foils, resulting in smaller 
lift, drag, and moment coefficients than for either of the 
smaller foils. It is interesting to note that this overall 
effect is contrary to blockage effects experienced in 
subcavitating water tunnel or wind tunnel testing, where 


force and moment coefficients are increased. 
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Baker also questions the sharp increase in wall effect 
predicted by Wu, Whitney, and Lin [6] for thin. symmetric 
wedges or for small attack angles in the lifting flow case. 
The results of the present work do not show a marked 
increase in wall effect on force and moment coefficients at 
the lower angles of attack, in agreement with Baker. 

Conversely, Baker's results indicate that the wall 
effect on the cavity plus wake can be tremendous. In the 
current research, this is, in fact, the most significant 
departure of the large foil data from the theoretical pre- 
dictions and from the data of the two smaller foils (Figs. 28 
through 30). 

Several attempts were made to bring the large foil data 
into closer agreement with the small and medium foil data by 
using existing two-dimensional blockage corrections. Wu, 
Whitney, and Brennen [4] derive wall effect corrections 
for the two-dimensional pure-drag (wedge) case, employing 
the open-wake and Riabouchinsky models. Wu, Whitney, and 
Lin [6] consider wall effects in two-dimensional lifting 
cavity flows, a situation more applicable to the current 
research; Baker [8] utilizes a computer program to translate 
this into information which could be readily applied to 


correct experimental data. 
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The pure-drag, open-wake model employs corrections to 
both drag coefficient and cavitation number. The drag 


coefficient corrected for blockage, C,', can be expressed as 


l+o' 2 
Cy’ = Cy tO 'CAe))., (1) 
1l+o 
c| 


where o', the cavitation number corrected for blockage, is 


given by 


leo | 2 
ok-= Ce = ‘ Cy) + OX), (2) 


Se 


where ope is the measured cavitation number, C, is the drag 


D 
coefficient as defined in the nomenclature, and i is the 
ratio of foil frontal width to tunnel width, which may be 


expressed as 
(mean chord) x sin(a) 


re (3) 
tunnel width 


For the large foil with A defined in this fashion, A = 0.30 
Sin (a). Application of the above corrections to the large 
foil data shifted the data points down and to the left on 

the plots of oy Js vs. C/O (Fags. 12 through 19) ‘but did 
not bring the data into any closer agreement with the small 


and medium foil data. When the cavity length data for the 
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large foil were replotted using the corrected cavitation 
numbers, the data points in Fig. 30 were moved downward; 
this downward movement was not sufficient, however, to bring 
the large foil cavity lengths into agreement with either 
theory or the small and medium foil data. 

Actually, the straightforward application of such a 
two-dimensional correction to the case of a finite-span 
wing should overcorrect for blockage, making \ much larger 
than it should be. Perhaps a better representation of i 


for the finite-span case would be 


foil frontal area 


3D tunnel cross-sectional area (without foil) 


(4) 
(foil planform area) x sin (a) 


tunnel cross-sectional area (without foil) 


Recalculating the corrections using 3p in place of A in 
equations (1) and (2) would make the corrections even smaller 
than previously calculated. Since the previous corrections 
proved to be inadequate, there is no reason to believe that 
the modified corrections should be any better. 

To correct measured cavitation number and measured 
drag coefficient utilizing the Riabouchinsky model, it is 
necessary to obtain Pp: the minimum pressure reading along 


the wall in the region of the cavity. Equation (1) is again 
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used to correct the drag coefficient, but now the corrected 


cavitation number is calculated as 


o' = (2/3)0,+ (17 3)G)s (5) 


where o'' is the cavitation number based on Por Pp and 
maximum velocity, V, and where oe is the cavitation number 
based on Por Pos and upstream velocity, U. It had orginally 
been hoped that a reasonably close measurement of the minimum 
pressure could be obtained by using only the downstream tap 
(tap D in Fig. 3); as was shown earlier, the interaction of 
the static pressure controller precluded making this measure- 
ment accurately. In an attempt to apply the correction based 
on the Riabouchinsky model, the gross assumption was made 
that the maximum frontal area of the cavity was approximately 
the same as the frontal area of the foil; utilizing a 
continuity argument and the steady-flow Bernoulli equation, 
the pressure minimum was estimated, thus allowing application 
of the wall effect correction based on the Riabouchinsky 
model. Unfortunately, this procedure gave results which 

were comparable to those obtained using the open-wake model. 
The failure of both models in the current application 

should not be surprising since they were designed for two- 
dimensional, pure-drag cavity flows, and not for three- 


dimensional, lifting cavity flows. 
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Whereas the pure-drag wall effect corrections changed 
both drag coefficient and cavitation number, the results 
of Baker's computerization of the two-dimensional Wu, Whitney, 
and Lin [6] lifting case wall effect model are presented 
in terms of a change to the force coefficients (i.e., 
drag and lift coefficients in this case), referenced to the 
uncorrected measured cavitation number, O,- Accordingly, 
Baker's corrections to the lift and drag coefficients are 
expressed as a percentage difference from the free stream 
condition. He presents data for only one ratio of tunnel 
height-to-arclength (i.e., chord length), H/c = 4. 
Calculating this ratio using the mean chord of the large foil 
gives H/c = 3.33; recalculating H/c as the ratio of tunnel 
cross-sectional area to foil area, H/c = 4.48 is obtained. 
It was concluded, therefore, that the H/c = 4 corrections 
given by Baker should be suitable for correction of the large 
foil data. Upon application of these additive corrections 
to both lift and drag data, it was found to bring the large 
foil data into close agreement with the small and medium 
foil data. 

Because an infinite wake model is utilized in 
reference [8], Baker's results give no correction for the 
extremely large cavity lengths recorded for the large foil. 


His results do show, however, that as blockage increases, 
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the nondimensional cavity thickness increases. This was 
borne out nicely when the cavitation characteristics of the 
three foils utilized in the current research were compared 

at an angle of attack of 8 degrees. On the small foil, 

which produced the least blockage, supercavitation was never 
totally achieved, no matter how low the static pressure was 
brought. For the medium and large foils, which produced more 
blockage than did the small foil, cavity termination was 


achieved behind the foil over a range of static pressures. 
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CONCLUSIONS 


As a general rule, the agreement of the experimental 
data with unbounded flow theory was good for the small and 
medium foils, whose ratios of half-span to tunnel width 
were 1/4 and 1/2, respectively, and whose ratios of tunnel 
height to mean chord were approximately 10 and 5, respectively. 
It appeared that in order to bring the small foil and medium 
foil data into agreement with theory, it was necessary to 
apply only standard wind tunnel procedures in correcting 
for the effects of the images of the trailing vortices. 

Of the three geometrically similar aspect ratio five foils, 
the largest foil, which protruded 3/4 of the way into the 
tunnel, and which had a ratio of tunnel height to mean chord 
of approximately 3.33, had force and moment data which were 
close to the small and medium foil data, although slightly 
lower; the cavity length data for the large foil differed 
greatly from the small and medium foil data, showing much 
longer cavities for the same ratio of cavitation number to 
angle of attack. It is evident that blockage corrections 
are required in addition to the standard wind tunnel down- 
wash corrections. 

Application and modification of blockage corrections 
based on various two-dimensional pure-drag cavity flow models 


failed to bring the large foil data into agreement with the 
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data for the two smaller foils. It was found, however, that 
corrections based on a two-dimensional lifting cavity flow 
model brought the large foil force coefficient data into 
much closer agreement with the small and medium foil 
data. While these force coefficient corrections were found 
to give an "engineering order of accuracy," no corrections 
have been found which give adequate corrections to the large 
foil cavity length data; it is clearly evident that further 
analytical work is necessary in this area. 

It is recommended that for future experiments, a static 
pressure measuring system similar to that depicted in Fig. 3 


be utilized. 
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(b) MEDIUM FOIL: @=12°, O,= 3473, Le =1.55 


(Cc) LARGE FOIL: @=12°, O,=.3265, L/c=21 
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APPENDIX I 
Standard Wind Tunnel Corrections for Effects of 
Images of Trailing Vortices 

The experimental data shown in Figs. 4 through 30 are 
corrected for the apparent upwash due to the interference 
of the tunnel walls, in accordance with standard wind 
tunnel practice. The boundary condition of zero normal 
velocity at the tunnel walls is satisfied in the analytical 
treatment of the problem by introducing a doubly infinite 
system of image vortices; the interference due to the 
presence of the tunnel walls is then calculated to be the 
interference caused by the series of vortex images on the 
flow at the model. The net result is an additive correction 
to the induced drag and the angle of attack (to be perfectly 
correct, there is also a slight negative correction to the 
measured lift, but this is usually neglected). Upon 
application of the correction factor, the flow about the 
model is transformed to the flow which would be seen if the 
foil were in an infinite stream. 

Since the hydrofoils tested in the square tunnel were 
only of half-span, a full-span model would require a duplex 
tunnel (that is, one whose ratio of height to width is 1/2). 
The ratio of model span to tunnel width for the small, medium, 
and large foils is 1/4, 1/2, and 3/4, respectively. Entering 


Fig. 6.34 of Pope and Harper [1], 6, the tunnel correction 
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factor for the influence of the images of the trailing 
vortex system was found to be 0.126 for the small, 0.104 for 


the medium, and 0.092 for the large foil. If the ratio of 


the foil area to tunnel cross-sectional area is S/So then 


the change in angle of attack is 
Cs: (radians) (I-1) 
The maximum value for Aa, encountered in the experiments was 


1.47 degrees. When Aa. is added to the measured geometrical 


angle of attack, a, the "true" angle of attack is obtained 


Gn = a + Aa, - VE=2) 


(In addition to the above correction, the data reduction 
program accounted for the dynamometer shaft twist, which 


never exceeded 0.06 degree.) The change in induced drag is 


Ac, = 6 hy 
ay fe) 
= ba, Cy Ci=5) 
so that the total drag coefficient is 
(a = ¢C + AC (I-4) 
e Duncor i 
where Cy is the uncorrected drag coefficient calculated 


uncor 
from the raw experimental data. 


=752 
APPENDIX II 


Leehey's Matched Asymptotic Solution for 
Supercavitating Hydrofoils of 
Large Aspect Ratio 


1 A synopsis of Leehey's theoretical results, as derived 


| in references [9] and [11], is given. 


| Leehey's linearized steady-flow theory, which utilizes 
| the method of matched asymptotic expansions, is valid to 
| first order in angle of attack, and to second order in the 
| reciprocal of the aspect ratio. Gravitational effects were 
| neglected in the derivation. The analysis shows that the 
lift, drag, and moment coefficients, and non-dimensional 
_ cavity length (measured from midchord) are given by the 


following relationships: 


| TO 2siny-l 
| Cy = = —— ’ (EBS) 
| siny (l+siny) A(1l+siny) 
ae 2siny-1 
Cy = ac, = —_—_————_— 4l - ——_+},, (II-2) 
sin y (l+siny) A(1+siny) 
-40 1-2siny+2sin2y 
Cu = ens, oo) - ——_}, (II-3) 
. 1 (1+siny) A siny (1+siny) 
and 
2 8 a 2 
L/c = sec’y - 1/2 - _ Ch r (II-4) 
TaA Cc 2D 


where C the two-dimensional lift coefficient may be 


é 
Lop 


expressed as 


=7iG— 


TO 
Cr, ee (II-5) 
2D siny (1+siny) 
and where 
| 20 
| yY = arctan |—|, (II-6) 
% 
and 
P,7P 
¢.= —— (ii=7)) 


C> (1/2) pu2— 


Experiments performed by Leehey and Stellinger [9] showed 
| 


that the above predictions for lift and drag coefficients 


and for cavity lengths were good; it was evident from the 


experimental results, however, that lifting surface corrections 


would be needed to improve the moment coefficient prediction. 


} 
} 
} 
} 
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APPENDIX III 


Numerical Lifting Surface Theory for Supercavitating 
Hydrofoils of Finite Span 


A brief description of the numerical lifting surface 
theory developed by Jiang and Leehey [12] for supercavitating 
hydrofoils of finite span will be given. 

The steady flow theory is linearized for application 
to thin wings at small attack angles. Discrete sources and 
vortices are utilized in the representation of the physical 


model, and the coupled integral equations are reduced to a 


set of simultaneous algebraic equations. The cavity closure 


condition is applied stripwise on the cavity in an iterative 
process to obtain the desired cavitation number over the 
supercavitated planform. 

This numerical theory should prove to be more accurate 
than the asymptotic theory of Leehey [11]; an implicit 
assumption in the latter theory is that the cavity length is 
less than the span, so that for small 0 /a (i.e.., long 
cavities), the numerical theory should indeed give better 
predictions. It should be noted, also, that for the 
O/ Om = 1.0 curves in Figs. 4, 5, 7, and 8 of reference |9], 
there is an ever-present difference between the experimental 
data and Leehey's asymptotic theoretical predictions for 
lift and drag coefficients. For long and short cavities, 


however, there is an even stronger disagreement between data 


re - — 
\s Ste Spey gas 
: ile - ctegia 
33 
- is hed | 
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and Leehey's predictions for moment coefficient (Figs. 9 
and 10 of reference [9]). The numerical lifting surface 
theory was developed in an attempt to overcome these 
difficulties. 

The foil-cavity surface is collapsed to become a 
region in the (x-z) plane of the foil. The relationships 
for the jump in the tangential and normal components of 
the perturbation velocity upon crossing the foil surface 


and cavity projection surface are: 


u(x,z, 0) ="ulx,z, 0) = -y(x,Z), (III-1) 


G(s, Z)- (III-2) 


v(x,Z,‘0) = ¥(x ;z,.0) 


where y and q are the vortex and source distributions, 


respectively. The boundary conditions are: 


dy (x,z 70) ; 
= “a re on the foil ‘J 
v(x.z, 0) = a 2 ie wetted surface (AE E= 3) 
x 
E on the cavity z 
nea o,/2 boundary eee 
& (2) 
= cavity closure 
ate reset = 0 condition, (III-5) 


X_ (Z) 


Where 2£(z) is the cavity length, measured from the leading 


edge as a function of spanwise position, and X, (2) is the 


— s 


= = 


opr! + ei 14565 1? Son wt ec iad 


"actin sor - 


yacsth fm mg 


mm 17d 


i it? ~» Sleog’ ak 
i= ; 


, , 


4 
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chordwise coordinate of the leading edge as a function of 
spanwise position. 

The discrete vortex and source method is utilized in 
the formulation of a series of simultaneous equations. 
Upon subdividing the projected surface into small elements, 
a discrete bound vortex, a trailing vortex, and a source 
are located within each of these elements. The bound 
vortex is situated on the quarter-chord line of each element, 
and induced velocities are calculated for all elements at 
their midspan, three-quarter chord positions. The 
concentrated source is taken to be a constant distribution 
across each element at its three-quarter chord position. 
The boundary condition of fixed cavity pressure is imposed 
on each element at its midspan, quarter-chord position. 
Utilizing these boundary conditions, the following equations 


are obtained: 


gq = 
ee J a. aa) als -5 Q ae if ; ee er »M1 (III-6) 
4n i,j 7 oT a Ne 
Y 1 q Oo 5 a joe ear (1 Ti=7 ) 
= ee) eee a Oe em 1.8 
a 74 oad | a ieee QO for k > Mi. 
G35 
} —i hs,. = 0, jane eee NG (LEI-8) 
i a J 
where Aiske and Stee are the v and u velocities at the (k,2) 


Wem 
nu? 6 an epbs Diiee 


) om 2S.3:90Tey uv S9e' gs 


=30= 


-th control point induced by the unit strength vortex and 


source at element (i,j). For the (i,j)-th element, eee is 


the element length, Ye the discrete vortex strength, and 


j 
qi the discrete source strength. Ml and N are the number 
of elements on the foil along the chord and span, respectively. 


M is the number of elements along the chord in the cavity. 
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APPENDIX IV 


Data Reduction Program and Data Printout 


Beginning on p. 82, a listing is given of the FORTRAN IV 
computer program utilized in the reduction of the experimental 
“data. 

The program output is also given, beginning on p. 92. 

| For each foil at each angle of attack, a complete printout 

‘of the experimental data is given in raw form, with corrections 
for instrument tares, in reduced form, in nondimensional form 
“without any corrections for wall effect, and finally in 


-nondimensional form with corrections for the effects of the 


images of the trailing vortex system. 
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